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VO Abstract In a recent study we demonstrated the emergence of turbulence in a trapped Bose-Einstein conden- 

sate of **^Rb atoms. An intriguing observation in such a system is the behavior of the turbulent cloud during 
free expansion. The aspect ratio of the cloud size does not change in the way one would expect for an ordinary 
C/3 non-rotating (vortex-free) condensate. Here we show that the anomalous expansion can be understood, at least 

^ qualitatively, in terms of the presence of vorticity distributed throughout the cloud, effectively counteracting 

^p* the usual reversal of the aspect ratio seen in free time-of-flight expansion of non-rotating condensates. 

Keywords atomic quantum fluid, condensate expansion, vortices, turbulence 

^ PACS 03.75.Hh • 05.30.Jp • 67.40.Db 

CT' 

g 1 Introduction 

Atomic quantum fluids have become an important system in the investigation of many effects inherent to 
P5 superfluidity UH^KBIT, 5 1 . One of the most striking properties of a superfluid is the behavior with respect to 

O rotations, which in these systems usually exist in the form of quantized vortices Ii6,lLi8ji9J . Vortices are among 

, ^ , the most investigated excitations in superfluids, and since their first appearance in atomic systems lll0lllll[T2l . 

they have also become an import aspect of research in atomic quantum fluids, 
psj Early experimental studies of a single vortex in a large trapped condensate, obtained by stirring techniques, 

^ allowed the first investigations of vortex nucleation and stability in these systems II13I . Subsequently, the 

On production of a doubly charged vortex state fl4], followed by the observation of its splitting into two singly 

charged vortices, has shown the dynamical instability of these multiply-charged configurations and confirmed 
that the most energetically favorable configuration contains only singly quantized vortices. As result, in the 
limit of large rotations and low dissipation, the vortices will equilibrate into an ordered lattice structure, 
similar to a crystal. 

Vortex lattices have been observed in atomic Bose-Einstein condensates (BEC) by several groups 1121 
I [T3]|15l[T6l . Those experiments have shown a broad range of vortex numbers A'v, from just a few to over a 
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hundred structures, almost always producing patterns similar to the Abrikosov lattices, observed in type-II 
superconductors. More recently, vortices have also been produced in trapped quantum gases by an excitation 
field that simultaneously combine motion and rotations along multiple axes [17, 18|. Numerical studies using 
the Gross-Pitaevskii equation (GP) have shown |19 20| that similar configurations evolve dynamically into 
the quantum turbulent regime, where tangled vortices are distributed randomly in the condensate. 

Many important aspects of turbulence in quantum fluids have been studied in the seminal work by Vinen 
1211 and others in superfluid helium. An important motivation to study quantum turbulence has been the un- 
derstanding of turbulence in a classical fluid, which is a very challenging and important practical problem. The 
unique restriction of motion of a superfluid, due to its quantum-mechanical constraints, results in a reduced 
number of degrees of freedom and makes the study of turbulence in quantum fluids (quantum turbulence) 
an attractive route to gain insight into classical turbulence 122) 123 Il24]| . Turbulence in atomic condensates, 
therefore, creates new exciting possibilities in this important field of research. 

In this communication, we present a simple hydrodynamical model that qualitatively explains the anoma- 
lous free expansion observed experimentally in [18|. Our model describes a condensate with vorticity dis- 
tributed along different directions, and uses the hydrodynamical description, in the Thomas-Fermi limit, to 
obtain analytical expressions for the time evolution of the aspect ratio of the expanding cloud of atoms, after 
its release from the trap. 



2 Production of a turbulent BEC 

In order to study a turbulent superfluid in the laboratory, we first produce a **^Rb condensate in an axisymmet- 
ric harmonic magnetic trap, with (Ox = (Oy = (Op = 2n x 207 Hz and (0, =2n x23 Hz, using the procedures 
described in reference 1251 . This produces typically a condensate with 2x10^ atoms which is then kept in an 
elongated (cylindrical) magnetic trap, where an external oscillatory magnetic field is applied, slightly off-axis 
and for variable amounts of time, by a pair of quadrupole coils. The resulting time-dependent field generates 
an oscillatory potential which imparts spatial and time varying forces to the atomic cloud. That combination 
produces simultaneous small rotations, translations and deformations of the condensate. 

Within a certain range of parameters, this perturbation leads to the formation of singly quantized vortices 
in the condensate, as the one present by Fig.[T^. Varying these experimental parameters one can produce very 




Fig. 1 Density profile of the atomic condensate cloud, obtained using absorption imaging after 15 ms of time-of-flight free 
expansion (images have different scales), showing different regimes: (a) single vortex in an ordinary condensate, where the 
Field-of-View (FOV) is 60 [im x 80 [im; (b) density distribution of a turbulent condensate, with FOV of 300 [im x 100 [im. 

complicated patterns, as seen in Fig. [TJ), which are reminiscent of the density distribution of a vortex tangle 
II 191 122 1. These density fluctuations remain in the cloud for times longer than several tents of milliseconds 
(longer than it takes for a typical density notch to move across the entire cloud), supporting the idea that these 
density fluctuations are due to tangled quantized vortices, and not simply a transient sound field caused by 
the excitation coils. Due to the limitations of the absorption imaging technique, which provides an integrated 
absorption along the line of sight, one cannot extract the full 3D arrangement of the vortices from images 
collected along a single direction. Comparison of the two density profiles in Fig.[T] however, gives a clear dis- 
tinction between the regular (a) and the turbulent (b) vortex regimes, based on the higher density fluctuations 
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present in the latter and the significant difference in the cloud sizes. But, as we discuss in the next section, an 
even more striking effect is the anomalous expansion observed ifTSl for a condensate in the turbulent regime. 
It is this behavior in time-of-flight that motivated the work presented here. 

The exact mechanism of vortex nucleation resulting in the turbulent regime is currently under study. We 
believe that it may be related to the scheme proposed in [19,20|, where the authors show that numerical 
simulation of the GP equation reveals turbulence induced in a trapped BEC by combined rotations around 
two axes (for more details about our trap excitation and the generation of the turbulent regime, see ref. Il26ll ). 



3 Observation of a self-similar expansion in a turbulent condensate 

One of the clear signatures of quantum degeneracy in trapped ultracold gases is the asymmetric time-of-flight 
expansion of the atoms released from an anisotropic trap (see, for example [6|, Chap. 7). This behavior is very 
distinct from a classical (thermal) gas, that always expands isotropically (unitary aspect ratio) at long enough 
time-of-flights. Quantum gases in anisotropic potentials will always have asymmetric velocity distributions, 
even for the ideal case where the asymmetry comes directly from the uncertainty principle. More often though, 
as it happens in most experiments involving alkali atoms with repulsive interactions, the effect is enhanced 
by the non-linear interactions, leading to an excess kinetic energy in the direction of tightest confinement. 
That happens because of the conversion of the mean-field energy along that direction. This behavior is well 
established in the field of quantum gases and has been studied and discussed by various authors [27 , 28 1, since 
the early days of BEC in alkali gases, being the general case typically well described by the GP equation. 

Therefore, for an initially (/ = 0) elongated vortex-free condensate with Thomas-Fermi radii Rp{0) ^ 
Rz{0), the expansion is much faster in the radial (tighter) direction. Consequently, the aspect ratio Rp (f ) /Rz(t) 
quickly grows in time-of-flight, from its initial small value and reaching an asymptotic value that depends 
only in the initial aspect ratio and the mean-field energy of the condensate (typically this final value is large). 
Hence, the free expansion of a non-rotating condensate produces a characteristic 'inversion' of the initial small 
aspect ratio, becoming larger (or smaller, depending on the directions) than unity at long time-of-flights. In 
contrast, the asymptotic value of the aspect ratio of a thermal (non-condensed) cloud always tends to the unity. 

A remarkable behavior of the turbulent quantum degenerate cloud in free expansion is that it seems to 
keep the cloud aspect ratio essentially unchanged. In fact, the measurements in ref. II18I show a complete 
suppression of the aspect ratio inversion within the range of time-of-flights observed. 

Here, we call this behavior self-similar expansion [29 30], and we explain it qualitatively in terms of a 
distribution of vortex lines, using a very simple model primarily intended to gain some intuition about the 
much more complicated turbulent regime. These three regimes are shown in Fig. |2] where one clearly sees 
that the turbulent cloud expansion (c) is markedly different from both a non-rotating (vortex-free) BEC (b) 
and a thermal cloud (a). 



4 Theory 



We start with the hydrodynamic equations ([T]l and Q that describe the macroscopic dynamics of the con- 
densate in the same way the GP equation. In these equations, n is the atomic number density, 1^ the velocity 
field, 

Vtrap the confining harmonic potential, and the interatomic interaction appears in the positive coupling 
constant g = Anfi^as/m (with as the s-wave scattering length and m the atomic mass). 
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These equations are equivalent to the continuity equation and Euler's equation of a perfect (inviscid) fluid. 

Typically for a superfluid, they reduce to represent a potential (irrotational) flow, since = 0. In that 

form, these equations are useful to calculate the collective oscillation in condensates [jSlEl^ where they have 
been applied very successfully ||3I||32||28]| . 

However, in the presence of vortices (singularities) the right-hand side of Eq. (j2| is not zero, which sig- 
nificantly alters the relevant physics. For a sufficiently large number of vortices, it is appropriate to write 
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Fig. 2 Evolution of the aspect ratio of tlie cloud during time-of-flight, for different initial conditions: (a) shows a classical 
thermal (non-condensed) cloud expansion, for which the value of aspect ratio tends to unity; (b) an anisotropic expansion of 
an ordinary BEC, without vortices, where there is an inversion of the initial aspect ratio; (c) expansion of a condensate in the 
turbulent regime, where the aspect ratio seems effectively fixed. 



V X This approximation of distributed vorticity has been used very successfully in Refs. II33II34I to 

describe the dynamics of a condensate containing a vortex lattice. Note that this step introduces a rotational 
component in the velocity field of the condensate, such that one can write 1^ = v/ + Vs, with v/ representing 

the irrotational component and vg = Q x Here Q is the macroscopic angular velocity of the sample, and it 
corresponds to the Feyrmian expression given in Eq. ([3]l, with ny being the areal vortex density {ny = Ny /A^), 
derived when the singular velocity fields of the vortices are distributed uniformly in the plane of rotation A±. 

^ h nh 

n = — ny = — ny (3) 
2m m 

The aim here is to describe the free expansion of a BEC containing angular momentum distributed along 
different directions. To first order we neglect the interactions among vortices, considering only the macro- 
scopic effect of the rotation in each direction. We also do not include any decay mechanism of the vortices (or 
turbulence) during time-of flight, because of the rapid decrease of the density after release from the trap ||6] 
13. Starting from the general description of rotational hydrodynamics, we focus our attention on two simple 
limiting cases, in order to develop some intuition about the system. In addition, since we are seeking to derive 
analytical expressions, here we will implicitly use the Thomas-Fermi (TP) approximation, where the kinetic 
energy is much smaller than the interaction energy. 

Consider the following ansatz for the density n and for the velocity field ~^ ll34ll : 



^(^,0 = \v{b,{t)x^+by{t)y^ + b,{t)z^) +n{t) X 7^, (5) 

Here = ln/mcOj^ gives the TF radii in terms of the chemical potential ji, and the peak (central) density riQ 
follows from Eq. m, with the normalization to the total number of atoms A^. Equations (4) and (6) basically 
express the TF approximation, known to represent well the system [6,7| under these circumstances, while (5) 
gives a general description of the velocity field in a rotating condensate. The coefficients bi{t) are temporary 
unknowns to be determined in what follows. 
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As a first example, we consider a vortex array aligned along the symmetry axis, that is, Q. parallel to z 
axis (see Fig.|3]l. Substituting Eq. (j4| and Eq. ([5]l in the continuity equation (TJ, we find the following relation 
between the general coefficient in the irrotational velocity ansatz bi{t) and the radius /?/(f) : 



Riit) 

Hty 



^'(0 = ^,;=w (7) 



Substituting Eq. (|7]) in the Euler equation (|2]), together with the Feynman expression Eq. ([8|, we find 
the evolution equations Q for the condensate radii (7?^ = Ry = Rp). For clarity, we simplified the notation 
(dropping the explicit time-dependence) in equation (|9jl, and used Rp to represent the radial (transversal) size. 



^2(0 = -^ (8) 
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(9) 



The numerical simulation of equation (j9| was done with experimental parameters of Ref. |18|, using a 
Runge-Kutta Ath order method. Starting from the equilibrium configuration (^,(0) = 0), we solved Eq. (nI 
without the trap term to obtain the evolution of the aspect ratio during the free expansion, as shown in Fig.p] 
for various values of initial vorticity £1. 
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Fig. 3 Representation of the system geometry (left) and evolution of the aspect ratio during time-of-flight (right), given by the 
Eq. (9), for different values of vorticity in a vortex lattice aligned along the symmetry axis z- 



It is clear that the extra kinetic energy due to the vorticity produces not only a larger initial aspect ratio, 
but also an extra acceleration of the radial expansion (perpendicular to £2). At higher values of Q this signif- 
icantly reduces the time necessary for the inversion of the initial aspect ratio. It is important to note that we 

— 1/2 

constrained the maximum value of Q used in our simulations, to maintain the vortex core size (^ «q ) 
small compared to the vortex line separation (y^h/niQ). Another way to verify the limiting value for £2 is 
comparing the effective energy due to the vorticity with the chemical potential: jU > HQ . The latter has to be 
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large in order to be consistent with the the TF approximation. Particularly, for the hmiting cases considered 
here, we have: 

Nv = ^ 11 = nns), nn = o 

Nv = 10^ 11 = 6h(o, hO. = Ihcb 
where both quantities, Q. and /i, were calculated at f = with the initial radius 7?, (0), and they are written in 

1 /3 

terms of the geometric average trap frequency 0) = ((Op^ft)^) 

The results of this first configuration illustrate clearly the contribution of the vorticity term, parallel to 
the symmetry vorticity. We now move on to consider the case where the vortex lattice is perpendicular to the 
symmetry axis of the trap {x-axis). This is, in fact, the preferred direction for the vortex lines in our experiment 
and therefore more representative of what we see in the laboratory. 

Again we derive equations for the evolution of the condensate radii, subject to this particular situation. 
Here, due to the asymmetry of the trap potential in the plane perpendicular to the axis of rotation, we included 
an additional irrotational term to the velocity field ansatz Eq. to produce an initially stable configura- 
tion 1351 . 

^C^.t) = i V {b,{t)^ + by{t)y^+b,{t)z^)+t{t)x^ + a{t)f{yz) (10) 



Similarly to the previous case, through the continuity equation we find the relations (|7|, ( [TT| and ( 12 1 for 

bi{t), Q{t) and a{t) respectively: 

X2(0 = ^— (11) 
mRy{t)R,{t)' 



Substituting into Eq. (|2|, we obtain the dynamical equations that describe the free evolution of the radii for 



this asynmietric configuration Eq. ( 1 3 i. The last term on the left-hand side of Eq.( 1 3 1 is related to the repulsive 



interaction force, whereas the term of the right represents the contribution of the vorticity. 
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The corresponding results for the evolution of the aspect ratio are presented in Fig. |4] which summarizes 
the most important result of this paper, as it reproduces qualitatively the non-inversion of the aspect ratio, 
provided the vorticity is above a certain value. Another point to observe from the simulations is that the vortex 
term alters significantly the dynamics of the expansion. The repulsive mean-field interactions act mostly to 
establish the initial equilibrium configuration and the expansion for short times (because they involve four 
inverse powers of the radii /?,). In contrast, the centrifugal (vortex) term with the conserved number of vortices 



Ny effectively has three inverse powers of Rj (the right-hand side of Eq. ( 13 1). It decays more slowly than 
the repulsive local-density contribution during the free expansion. Hence the centrifugal term dominates the 
expansion dynamics at long times. 
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Fig. 4 Evolution of the aspect ratio during time-of-fligiit, for a system with vorticity along the the x-axis (perpendicular to the 
symmetry axis). The non-inversion of the aspect ratio (the radial divided by the axial dimension) is shown by the solid line of the 
left hand figure. The figure on the right shows the acceleration of the expansion in the plane perpendicular to the vortex array. 



5 Conclusions 

In this paper we explain qualitatively the anomalous expansion of a condensate in the turbulent regime. Our 
model shows that he behavior of the expansion depends intrinsically on the amount and distribution of vor- 
ticity present in the sample. To obtain this result we used a simple macroscopic description, based on the 
rotational hydrodynamics, of a condensate with a relatively large vorticity aligned both parallel and perpen- 
dicular to the symmetry axis of the trap. We believe that due to the geometry of the exciting oscillating field 
in our apparatus 1 18], more angular momentum is coupled into directions perpendicular to the symmetry axis 
of the trap. That may explain, at least qualitatively, why we obtain such a significant reduction of the aspect 
ratio inversion in the experiments. 

In order to describe quantitatively the results observed in the laboratory, it is clear that a more elabo- 
rate numerical simulation of the turbulent regime, accounting for the distributed vorticity and including the 
vortex dynamics and possibly reconnection effects (mostly in the trap, since this effect should be small in 
time-of-flight, when the density decays very rapidly), will be needed to fully account for the experimental ob- 
servations. The goal here was to present a simple model that seems to capture the main physical mechanism 
behind the intriguing anomalous free expansion of the cloud. Although very phenomenological in nature, this 
simple toy-model already gives a clear indication that a disordered volumetric distribution of vortex lines can 
greatly affect the dynamics of the density profile in time-of-flight. 

Acknowledgements We thank M. Tsubota for helpful discussions and E. A. L. Henn, J. A. Seman and D. V. Magalhaes for 
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